WUB/12-14 



Mean-Field Gauge Interactions in Five Dimensions II. 

The Orbifold. 



Nikos Irges 1 , Francesco Knechtli 2 and Kyoko Yoneyama 2 

1. Department of Physics 
National Technical University of Athens 
Zografou Campus, GR-15780 Athens Greece 

2. Department of Physics, Bergische Universitat Wuppertal 
Gaussstr. 20, D-^2119 Wuppertal, Germany 

e-mail: irges@mail.ntua.gr, knechtli@physik.uni-wuppertal.de, 
yoneyama@physik.uni-wuppertal.de 



Abstract 

We study Gauge-Higgs Unification in five dimensions on the lattice by means of the mean- 
field expansion. We formulate it for the case of an SU{2) pure gauge theory and orbifold 
boundary conditions along the extra dimension, which explicitly break the gauge symme- 
try to U(l) on the boundaries. Our main result is that the gauge boson mass computed 
from the static potential along four-dimensional hyperplanes is nonzero implying spon- 
taneous symmetry breaking. This observation supports earlier data from Monte Carlo 
simulations |12|. 



1 



1 Introduction 



The phase diagram of five-dimensional gauge theories is surprisingly rich. On an infinite, 
hypercubic, anisotropic lattice it is parametrized by the two dimensionless parameters 
/?4 = /3/7 and /3 5 = (3^, where (3 is the lattice coupling and 7 the anisotropy parameter. 
In part I of this work [lj we explored the phase diagram of a five-dimensional SU (2) gauge 
theory with fully periodic boundary conditions, using an expansion around a mean-field 
background [2]. We concentrated on the regime where /3 ~ 0(1) and 7 < 1 where a line 
of second order phase transitions was observed. In the vicinity of this phase transition the 
system reduces dimensionally to four dimensions and in the continuum limit the physics 
is consistent with what one would expect from the lightest states with four dimensional 
quantum numbers. For the fully periodic system this would be a four dimensional SU(2) 
gauge theory coupled to an adjoint scalar in the confined phase [3]. 

In this work we extend the construction of [1] by changing the boundary conditions in 
the fifth dimension from periodic to orbifold. The embedding of the orbifold projection in 
the geometry introduces boundaries at the "ends" of the fifth dimension (which is now an 
interval) and its embedding into the gauge group as well known by now [I], alters the field 
content surviving on the boundaries. In this respect, for SU(N) there are at least two 
possibilities. One is when the orbifold action is such that the adjoint set of scalars (i.e. 
the extra-dimensional components of the gauge field) is projected out at the boundaries 
and one is left there with just a pure SU(N) gauge theory. This construction allows one 
to carry out an analysis similar to [3] but in a context that is directly generalizable to 
QCD once fermions are also added. The other possibility is to use the orbifold action to 
project out some of the gauge fields and some of the scalars. With such a choice it is 
possible to realize a field content similar to the one of the bosonic sector of the Standard 
Model. At a first stage we consider an SU(2) bulk gauge group. The orbifold then leaves 
a f/(l) theory coupled to a complex scalar on the boundaries. This setup could serve as 
the simplest prototype of the Higgs mechanism. 

The idea that the Standard Model Higgs particle may be the remnant of an extra 
dimensional gauge field is not new |5J. Also, investigations of five- dimensional gauge 
theories with a lattice regularization have both an analytical and a Monte Carlo past. 
The analytical work has been concentrated around the question of the existence or not 
of a layered phase [6] and around the existence or not of an ultraviolet fixed point [7]. 
Monte Carlo simulations have been mostly looking for the layered phase [8] , for first order 
bulk [9] (the pioneering work in this direction) or second order [10J (finite temperature 
and bulk) phase transitions and dimensional reduction via localization [11] . All of these 
lattice investigations have been carried out on fully periodic lattices. Recently there has 
been interest also in lattices with orbifold boundary conditions [12], [13]. Our motivation 
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to look more carefully at the phase diagram of five- dimensional orbifold gauge theories 
stems from the fact that Coleman- Weinberg computations [H], continuum perturbation 
theory at one loop [15] and exploratory lattice Monte Carlo simulations p2] indicate that 
once a Higgs mass is generated by quantum effects, it seems to remain finite, despite the 
non-renormalizable nature of the higher dimensional theory. 

A first attempt to probe analytically the regime away from the perturbative point, in 
order to see if there is a dynamical mechanism of spontaneous symmetry breaking (SSB) 
triggered by the gauge field (without the presence of fermions) was made in [16]. The 
idea there was to use the SU(N) Symanzik lattice effective action [17J (a is the lattice 
spacing and N 5 the number of lattice points in the fifth dimension) 



which is determined by a finite number of dimensionless coefficient functions c^ Pi \N^, (3) 
on an infinite spatial isotropic lattice, provided that one can consistently truncate the 
expansion. The ansatz in [16] was to truncate the expansion after the first two higher 
dimensional operators: one of order a corresponding to the lowest order dimension five 
boundary counterterm parametrized by the coefficient function c^(N^, (3) and one of 
order a 2 corresponding to the dimension six bulk operator tr(-DF) 2 parametrized by the 
coefficient function (N 5 , (3). To extract information about SSB in this setup, one 
assumes a vacuum expectation value (vev) v for one of the A 5 components of the five- 
dimensional gauge field Am and uses the truncated expansion expanded around this vev 
to compute a Coleman- Weinberg type potential V{pi) for the dimensionless quantity 



with g 5 the five- dimensional coupling, ttR = N 5 aa the size of the fifth dimension and 
we have written this formula for an anisotropic lattice with lattice spacings 0,4 and 05 
(7 = 04/^5 in the classical limit). One then finds the preferred value for a by minimizing 
V(a) rjUdH] and calls it a m i n . As a consequence, the otherwise massless gauge boson 
develops a mass due to this vev equal to 



This is the Hosotani mechanism, applied to the case of the orbifold. The result of the 
analysis of [16], performed at 7 = 1, was that indeed there exist values of and that 
yield a "Mexican hat" Higgs potential that triggers SSB with the Higgs particle having a 
mass of similar order as the gauge boson. In particular, it was shown that a non-zero 
is able to trigger SSB by itself, by shifting a min from integer (for which there is no SSB) 




(1.1) 




(1.2) 
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to half integer. This could be the main phenomenological gain from the complications 
encountered by entering in the interior of the phase diagram, in view of the fact that in 
the conventional continuum approaches where one takes c^\N 5 ,[3) = 0, one necessarily 
needs fermions in order to trigger SSB p2] (i.e. a non-integer a m i n ) and even if SSB is 
achieved, the Higgs typically turns out to be generically too light [19]. In the absence 
of a non-perturbative control of the theory, in [16] the coefficients c^ 5,6 ^ were treated as 
free parameters. Non-perturbatively however they are not free parameters and it is not 
guaranteed that the quantum theory generates values for the coefficients that trigger SSB. 

In this work we improve on these approximations by computing the Wilson loop and 
the mass spectrum of the lightest states, in the mean-field expansion far from the five- 
dimensional perturbative point, close to the bulk phase transition. The formalism involved 
is very similar to the one developed in [1] and therefore will be heavily used. We show 
that the mean-field expansion predicts the spontaneous breaking of the boundary gauge 
symmetry already in the pure gauge system and allows for a Higgs like scalar of similar 
mass as the mass of the broken gauge field. With infinite four- dimensional lattices, the 
parameters of our model are /3, 7 and N$. To the extent that the mean- field expansion is a 
good description of the non-perturbative system, any result stemming from this approach 
should be taken seriously. In fact, the first exploratory Monte Carlo studies of the orbifold 
theory [12] had earlier reached similar conclusions. 

In Section [2] we give a short review of the mean-field expansion formalism. In Section [3] 
we apply the general formalism to the five-dimensional SU(2) lattice gauge theory with 
orbifold boundary conditions along the extra dimension. In Section H] we present our 
numerical results and in Section [5] our conclusions. In the Appendices we detail the 
mean-field calculations of the propagator with orbifold boundary conditions and of the 
mass spectrum. 

2 A short review of the mean-field formalism 

The partition function of a gauge theory on the lattice is 



where <SV[C/] is the Wilson plaquette action and p denotes oriented plaquettes (i.e. each 
plaquette is counted with two orientations). In the mean-field approach [2] the link 
variables U are traded for the complex quantities H and V and in terms of these one 
rewrites the partition function as 



Z= BV DHe' s ^ [v ' H] , S eS = S W [V] + u{H) + {l/N)Retr{HV} , (2.2) 




p 



(2.1) 
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where the effective mean- field action u(H) is defined via 



(2.3) 



The mean-field or zeroth order approximation amounts to finding the minimum of the 
effective action when 



H 



HI, V — >V1, S eS [V, H] =minimal . 



(2.4) 



The zeroth order saddle point solution or "mean-field background" [V] can be easily 
obtained by taking derivatives of Eq. ( 12. 2 p with respect to V and H and require them to 
vanish. One then has 



v ■ () " 



OH 



7I= dS w [V] 



H 



dV 



(2.5) 



The above two equations are the ones that make the action extremal and define the 
mean-field solution to zeroth order. The free energy per lattice site is 



F=-^l„(Z). 



At O'th order we simply have 



S cS [V,H] 



Gaussian fluctuations are defined by setting 

H = H + h and V = V + v. 



(2.6) 



(2.7) 



(2i 



We impose a covariant gauge fixing on v. In [2U] it was shown that this is equivalent to 
gauge-fix the original links U. The integral 

(27r)l ft l/ 2 (27r)H/2 



y/det[(-l + PWj] 



introduces the pieces of the propagator 



5 2 S, 



cff 



5H 2 

'cff 



V,H 



h 2 = fnK^hi = h T K^h 

1 v J 



s 2 s c 



5V5H 
5 2 S eS 



V,H 



vh = v t K\f ) h 1 =v T K^h 



SV 2 



2 ts( vv ) T ts(vv) 

v = ViK-j 'vj = v K K >v 



(2.9) 

(2.10) 
(2.11) 
(2.12) 



which will be used extensively later. The quadratic part of the effective action is [v,h] = 

1 {hT K {hh) h + 2v T K (vh) h 

+ v T K^ vv ^v) , \h\ and \v\ denote the dimensionalities of the fluc- 
tuation variables h and v. 
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We would like to compute the expectation value of observables 

(O) = ^ J BU 0[U]e- Sw[u] (2.13) 
in the mean-field expansion. To first order it is given by the formal expression pQ 



1 I ft 2 CI 

with 



K~ x } , (2.14) 

v 



K = -K {vh) K {hh) 1 K {vh) + KW + if (gf) (2.15) 

and the second derivative of the observable is taken in the mean-field background. is 
the contribution from the gauge fixing term. actually turns out to be proportional 

to the unit matrix and drops out from all expressions. The free energy at this order 
becomes 

F (i) = j p(o)_i. ln ^_ (2.16) 

To extract the mass spectrum, we denote a generic, gauge invariant, time dependent 
observable as 0{t) and its connected version as O c {t) = O(t + t)O(t ). Defining the 
correlator 

C(t) =< O c (t) > - < O{t + 1) >< O(t ) > , (2.17) 
to first order in the fluctuations the expression reduces to C(t) = C^(t) with 

cm{t)= Kj^m K -A : „ 18) 



2 { 5 2 V 

where the notation t^ 1 ' 1 ) means one derivative acting on each of the O(t + 1) and O(to). 
The mass of the lowest lying state is then 

C (1) (t) 

m = lim In — . (2.19) 

It turns out that in order to extract the mass of the vector one needs to go to second 
order in the mean-field expansion. Physical expectation values are formally given at this 
order by [1] 

+ iwi) ((K^K-^ + iK-^iK-^ + iK-^K- 1 )^. 

i,j,l,m ^ ' V lm 

(2.20) 
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To extract the mass from the connected correlator is straightforward. Again, all time 
independent contribution (self energies) cancel from connected correlators and at the end 
the mass is obtained from 

m limln g"'W +<**(') (22 1) 
where C^'it) is the next to leading order correction to the Z-boson correlator. 



3 The lattice orbifold in the mean-field expansion 



We will now apply the formalism we described in general terms to a specific example: an 
SU{2) lattice gauge theory in 5 dimensions with Dirichlet boundary conditions for certain 
components of the gauge field along the fifth dimension. 

The discretized version of the S l /Z2 orbifold defined on a five- dimensional Euclidean 
lattice was constructed in [21]. The points on the lattice are labeled by integer coordinates 
n = {tim} but we will often use the notation Uq = t for the time component. The periodic 
spatial directions (M = \x = 1, 2, 3) have dimensionless extent L = Z/a 4 and the time-like 
direction (M = 0) has extent T. The fifth dimension (M = 5) has extent N$ = TcR/a^. 
The gauge-unfixed anisotropic mean-field Wilson plaquette action reads 



«M "5 = 1 

4ee 



N5-1 



cff 



n u "5=0 



^ tr V p Abound {n; fl, v) 

p,<v 

^ Re tr V^bound(n; fi, 5) 



4e 



^ Re tr V^gbound (n; fx, v) 

,H<v ri5=0, N$ 



JV5-I 

+EEE 

rifi ",5=1 /j 
JV5-I 

+EE 

rip "5=0 



+EE E 

«M M n 5 =0,iV 5 



u 2 (p(n, //)) + ^ M n > / i ) u a( n > Z 1 ) 
u 2 (/£?(n, 5)) + /i a (n, 5)u a (n, 5) 

«i(p(n, n)) + ^ h »( n ' l^) v a(n, n) 



(3.1) 



where the effective mean-field actions u\ and 112 are computed in Appendix |A] The cou- 
plings are defined as 

A= ^, &= ^. (3.2) 



2 ' 

si 



#505 
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In this work we will parametrize the infinite anisotropic lattice by the parameters f3 and 
7, where (3^ = and j3 5 = ft'y. A gauge transformation acts on a bulk link as 

U(n, M) — ► rt su ^ (n)U(n, M)OS su ™ (n + M) (3.3) 
on a boundary link as 

U(n, M) — ► (n)U(n, M)fi (l/(1))t (n + M) (3.4) 

and on a link whose one end is in the bulk and the other touches the boundary as 

U(n, M) — -> fi (C7(1)) (n)U(n, M)fi (5C/(2))t (n + M) . (3.5) 

One possibility is to derive the orbifold theory from its parent circle theory. In this setup 
a general link satisfies the orbifold projection condition 

T U(n, M) = U(n, M), r = T g K (3.6) 

where the reflection property about the origin of the fifth dimension is 

TZU(n,fi) = U(n,fi) 

KU(n,5) = U\n-5,5) (3.7) 

with 

n = (n jU ,n 5 ), n = -n 5 ). (3.8) 
The transformation property under group conjugation is 

T g U(n,M)=gU(n,M)g' 1 . (3.9) 

The boundary conditions that the above projections imply at their fixed points amount to 
Dirichlet boundary conditions for some of the links at the orbifold boundary hyperplanes. 
Consequently the gauge group variables at the boundaries are restricted to the subgroup 
of SU(N) invariant under group conjugation by a constant SU(N) matrix g with the 
property that g 2 is an element of the center of SU(N). Only gauge transformations 
that commute with g are still a symmetry at the boundaries and thus there, the orbifold 
breaks explicitly the gauge group. In general, a bulk group G breaks by g to an equal 
rank subgroup % on the two boundaries @]. 

For SU(2), which is the gauge group of our focus, we will take g = — icr 3 . This 
means that at the orbifold fixed points the gauge group SU(2) is broken to the U(l) 
subgroup parametrized by exp(z0a 3 ), where <fi are compact phases. In the continuum 
limit this implies that A 3 ^ (the 11 Z gauge boson" ) and A 1 ^ 2 (the "Higgs" ) satisfy Neumann 
boundary conditions and A 1 ^ 2 and A\ Dirichlet ones. 
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In the mean-field approach, we parametrize the fluctuating fields in the bulk as 

3 

V(m,M) = v (n,M) +i^2v A (n,M)a A , 

A=l 

3 

H(m, M) = h {n,M) -iJ2 h A{n,M)a A . (3.10) 

A=l 

The a A are the Pauli matrices. On the boundaries instead, we use the parameterizations 

V(n,M) = Vo(n, M) + iv 3 (n, M)a 3 , 

H(n,M) = h (n, M) — ih 3 (n, M)a 3 . (3.11) 
with vq,a £ C For later convenience we define the line 

l (n5 \t ,m) = J] V((t Q ,m,m 5 )-5) (3.12) 

m5=0 

and introduce the matrices 

a a = {1, ia A }, W a = {l,-ia A }, A = 1,2, 3. (3.13) 

For the computation of the Z and Higgs masses we first define the orbifold projected 
Polyakov loop 

P<°>(t,m) =l {N5 \t } m)gl^{t,m)g\ (3.14) 

satisfying TP^ = p(°\ in terms of which we define the field <&(°'(t, m) = p(°\t,m) — 
p(°^(t,m) and then the displaced Polyakov loop [lj 

zf )A (t, m) = a A V{{t, m, 0); k) $ (0)t (t, m + k) V((t, m, 0); fc) f $ {0) (t, m) , (3.15) 

which assigns a vector and a gauge index to the observable appropriate to a gauge boson. 
The Higgs observable is derived from the averaged over space and time location connected 
correlator 

O c H (t) = E HP i0 \t ,m')}tr{pV\t + t,m")} (3.16) 

to rh',fh" 

and the Z-boson from the correlator 

°z(t) = ^EE EE^° ),j4 ^^ , )}tr{4 0),A (*o + ^'')}. (3.17) 

to ml ,rn" A k 

Out of the above defined objects one can straightforwardly extract the masses using the 
general results of the previous section. 
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3.1 The mean-field background 



In order to determine the background we need the effective potentials U\ and u^. The 
effective potentials in the bulk are the same as on the torus, while on the boundaries they 
are 



Ul {H{n,M)) = -ln{I {p)), p = y/(Re h (n, M)) 2 + (Re h 3 (n, M)) 2 . (3.18) 

For details see Appendix lAl 

Translation invariance along the dimensions \i — 0, 1, 2, 3 means that we can parametrize 
the saddle point solution, which minimizes S e g, as follows [22J: for n 5 = 0,1, . . . , N 5 (four- 
dimensional links) 

H(n,fi) = h (n 5 )l , V(n, y) = v (n 5 )l , Vn M ,/i, (3.19) 
and for = 0, 1, . . . , JV5 — 1 (extra-dimensional links) 

H(n,5) = h {n h + 1/2)1, V{n, 5) = v {n 5 + 1/2)1 , Vn M . (3.20) 
The action at zeroth order reads (Af = L 3 TAf$, A/5 = N$ + 1) 



S cS [V,H] 



'JVb-1 



£ + k(o) 4 + k(iv 5 ) / 



n 5 = l 



iVfi-1 



- 1) J2 ^oK)(^oK + l/2)) 2 v (n 5 + 1) 

n 5 =0 

M>-1 



+ (d-l) 

iV B -l 



n 5 =0 



+«i(/io(A r 5 )) + £ ^2(^0(^5)) + £ M^Hfas) 

"5 = 1 

+ £ M^oK + V2)) + ^oK + l/2)U (n 5 + 1/2)] 



(3.21) 



n 5 =0 



The minimization equations lead to the following relations: for 715 = 

*„(0) = -^(0)) = £§M, 

MM )) 

^o(O) = /3 4 [(^-2)(l;o(0)) 3 + 7 2 (^(l/2))%o(l)] • 



(3.22) 
(3.23) 



A prime on u\ or «2 denotes differentiation with respect to its argument. Similarly, for 
n 5 — A5 we have 



^0(^5 



-u[(h (N 5 )) 



h(h (N 5 )) 



(3.24) 



I (h (N 5 )) 

h (N 5 ) = /34{d-2)(v {N 5 )) 3 + j 2 v (N 5 -l)(v (N 5 -l/2)f] . (3.25) 
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For 7i5 = 1, . . . , N5 — 1 (four-dimensional links) 

v {n 5 ) = -u 2 {h {n 5 )) = , (3.26) 

h{h {n 5 )) 



h (n 5 ) = (3, 



2(d - 2)(U (n 5 )) 3 + 7 2 (iv (ns + l/2)) 2 v (n 5 + 1) 
+UoK-l)(^oK-l/2)) 2 )] . (3.27) 
For 7i5 = 0, . . . , N$ — 1 (extra-dimensional links) 

I7 (n 5 + l/2) = -^ 2 (Vn 5 + l/2)) = ^^ oK + 1 ^ , (3.28) 

h{ho{n 5 + 1/2)) 

h (n 5 + l/2) = 2(3 5 (d-l)v (n 5 )v (n 5 + l/2)v (n 5 + l). (3.29) 

3.2 Observables from fluctuations around the background 

In sect. 2 we described the general formalism for computing observables from fluctuations 
around the mean-field background. Here we apply this formalism to our case and give 
the results for the free energy, the scalar and vector masses. 

The lattice propagator, Fourier transformed along the spatial and time directions is 
an object that contains the information about the boundary conditions. We denote its 
components as 

R- 1 = K- l {p\n^M\a'-p",nlM",a") . (3.30) 

The momenta p',p" are four dimensional momenta, however we will further split the 
momenta into their time and spatial components: p = po,Pk whenever it is necessary, for 
clarity. For a more detailed computation of the propagator we divert the reader at this 
point to Appendix B. 



3.2.1 The free energy 



The free energy to first order is 

F m = F (o) + ±_ ln 

2Af 



II det (T(p', n' 5 , M', a; p", n'>, M", a)D 



a=0 



(orb) -2 
FP 



(3.31) 



where is the determinant of the Faddeev- Popov matrix, also computed in Appendix 

B. The matrix T is defined as 



(3.32) 



There are torons both in the matrix T and D^p b \ which can be regularized as on the 
torus pQ , see the end of Appendix [B] 
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3.2.2 The Higgs and Z -boson masses 



For the Higgs, there is a non-trivial contribution already at first order. The result is 

where is the Polyakov loop Eq. ( I3.14p evaluated on the background and (0, 0) is 
defined in Eq. (1C.2j) . The above correlator does not contain torons since the 00 component 
of the propagator does not contain any. 



For the Z the result is 



p 1 k 



where 11^^(1,1) is defined in Eq. flC2[) . It contains regularizable torons, that is simul- 
taneous zero modes in the propagator and the observable, whose contribution vanish in 
the infinite lattice volume limit. 

Derivations can be found in Appendix ICl 



3.2.3 The static potential 

On the orbifold we have the three types of static potentials. Here we will be interested 
in the potentials extracted from Wilson loops in the four- dimensional hyperplanes, along 
either one of the boundaries and in the middle of the orbifold (i.e. at n§ = N§/2). We 
consider the Wilson loops of size r along one of the three spatial dimensions and we average 
over the possible orientations. The exchange contribution (to 5 2 O c /SV 2 ) at n 5 = is 

O ex = ^2(^ (0)) 2 ( i+ ^- 2 Wm-o 

\M=1,2,3 / k=l 



(3.35) 



and the self energy contributions 



O sc = ^2(v (0)) 2{t+n3) ~ 2 hi>oSM>>o 



5n 5 o(S a ' 5 a "o - S a ' 3 5 a "3)5p' o o5 p ^ I J^J <V m pm I 2 ^o^ngo • (3.36) 

\M=1,2,3 / 
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As for the torus [i], to first order we would like to compute 



c w - 2 ^2 ^2 

a', a" p' k n' B ,n'^ 



O^pln^O^'^^lnlO^'^K-Ho^ln^O^'^^lnlO,^') , (3.37) 



where O = C ex + O se , which is to be substituted in the general expression for the first 
order corrected static potential 

1 

V = const. - lim %r . (3.38) 

t 0[V] 

Applied to the gauge boson exchange between two static charges on the boundary the 
general formula reduces to 

K(0) = -.o g( , (0)V^^E 
i [2 cos (p' k r) + 2] K- 1 (0,p' k , 0, 0, 0; 0,p' k , 0, 0, 0) 

+ ±J2[ 2 cos (^ r ) - 2 ] K ~ X (°> Pk,°A3;0, p' k , 0,0,3) I. (3.39) 

The formula for the static potential Va(Ns/2) for the static potential along the four 
dimensional hyperplane in the middle of the orbifold is similar to Eq. (13. 39 p . the difference 
being that the background and the propagator are evaluated at n 5 = iVs and in the last 
line of Eq. (I3.39P there is a sum over all the gauge components a = 1, 2, 3. 



4 Spontaneous symmetry breaking 
4.1 The phase diagram 

Using the equations that determine the mean-field background in Section 13. 1\ one can 
extract the leading order approximation to the phase diagram. The equations are solved 
iteratively and numerically. The confined phase is defined as the phase where ^0(^5) = 
= vo(n 5 + 1/2) for all n 5 . When vo(n 5 ) 7^ and vo(n 5 + 1/2) = for all n.5, we define 
the layered phase. The Coulomb phase is defined where vo(n 5 ) ^ and vo(n 5 + 1/2) 7^ 
for all 715. We do not find a phase where ^0(^5) = and ^0(77.5 + 1/2) 7^ for all n 5 . The 
background is sensitive to /3, 7 and N 5 . On Fig. [T]we plot the phase diagram with color 
code, red for the confined phase, blue for the layered phase and white for the Coulomb 
phase. Green is used where for some reason the iterative process does not converge to a 
solution. For the rest of this paper we will stay in the Coulomb phase. 
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Figure 1: The mean-field phase diagram of the SU(2) orbifold theory in the (/3,7, N$) space. 
The color code is explained in the text. 

As for the torus, the line that separates the Coulomb from the confined phase is of first 
order for 7 larger than a value which is close to O.7., below which it turns into a second 
order phase transition. The order of the phase transition that the mean-field predicts 
must be taken with care though, as a more careful, fully non-perturbative analysis should 
be done. 

Next, using the quantities , and C$ we analyze the physical properties of the 
system on this phase diagram, with an emphasis on the issue of spontaneous symmetry 
breaking (SSB). 



4.2 The Higgs 

The Higgs mass in units of the lattice spacing Mh = a^niH, extracted from Cg in 
Eq. ( 13 .331) . depends on /3,7 and N5. The physical quantity of our interest is the Higgs 
mass in units of the radius of the fifth dimension 

F 1 =m H R = M H — . (4.1) 

7 7T 

In perturbation theory, the one-loop result [15] for SU(N), expressed in lattice parameters 
(relevant for the isotropic lattice) is (C 2 (N) = {N 2 - 1)/(2N)) 

^r = "^. «=^C(3)ft(JV). (4.2) 

On the left plot in Fig. [2] we show the iVs-dependence of Mh for 7 = 1 at /3 — 1.677 near 
the phase transition. We can see clearly that the perturbative formula is not valid at a 
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1/N 



Figure 2: The Higgs mass Mjj (left plot) and the Z boson mass Mz (right plot) as a function 
of I/-/V5 at 7 = 1 for f3 = 1.677. The squares are mean-field data, the line on the right plot is a 
quadratic fit and the line on the left plot a linear fit. 



generic point on the phase diagram. The line on the left plot in Fig. [2] is a quadratic 
fit. The phase transition is of first order, which means that the mass in lattice units M# 
cannot be lowered to zero but approaches a non-zero minimal value, which at = 1.677 
is approximately 0.69. 



4.3 The Z boson 



The Wilson loop can decide if there is SSB. We will choose our lattices so that the system 
is dimensionally reduced to four dimensions. Then we can describe the boundary gauge 
theory in four-dimensional terms. If the boundary £7(1) symmetry is spontaneously broken 
then the corresponding static potential extracted from Cyy Eq. ( 13.370 should be fitted by 
a (Ad) Yukawa form rather than by a Coulomb form. Starting from 

e -m z r 

VAr) = -b + const. , (4.3) 

r 

where b is a constant, we define the quantity y(r) = log(r 2 F4(r)) where F±(r) = d V±(r)/dr 
from which we form the combination 

a 4 y'(r) = -M z + — ^- . (4.4) 
rrizr + 1 

By Mz = a^mz we denote the Z mass in lattice units. We then determine Mz iteratively 
by requiring that a plateau for —a^y'(f) + Mz/(rrizr + 1) forms. The plateaus, for large 
enough L, stabilize as L is further increased, so that Mz at infinite L depends on (3, 7 
and N 5 . Note that in the case of a first order phase transition this corresponds to a Z 
boson in an infinite physical volume at a finite lattice spacing. 
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Figure 3: The ratio of the Higgs to the Z boson mass Eq. (|4.5p . Comparison of Monte Carlo 
(diamonds [12] and circles [16]) and mean- field data (squares) at 7 = 1. 



4-3.1 Isotropic lattices 

On the right plot of Fig. [2] we show the Mz plateau values as a function of I/-/V5 at 
fixed = 1.677 and 7 = 1, near the bulk phase transition. The plateau values of Mz 
do not depend on L for L > 200 and there is no sign of a plateau for a zero mass. A 
linear fit with slope 3.32, which is very close to ir, describes the data very well. The 
most striking observation is that the boundary gauge boson is massive, pointing to the 
dynamical spontaneous breaking of the U(l) symmetry. Clearly, since (3 and 7 are kept 
fixed, the masses on Fig. [2] correspond to different lattice spacings (the location of the 
phase transition /3 C depends on N$). Nevertheless, from the Kaluza-Klein description 
we expect to see an approximate l/N 5 dependence. Comparing the data of Fig. [2] to 
Eq. (II. 3p our result indicates a value a min = 1. However the regime of Higgs masses 
in Fig. [2] corresponds to values of F\ = tuh R 3> 1 which is not the regime where the 
Coleman- Weinberg calculation is performed and for which Fi < 1. In perturbation theory 
the minimum at a min = 1 is equivalent to the one at a min = and describes a situation 
without SSB. Clearly this is not the case for the mean-field data. 

On Fig. [3] we plot mean-field data (squares) for the ratio 



Phz 



m z 



(4.5) 



obtained using N$ = 4,6,8 and L = 200. As far as we have checked the mean-field results 
in Fig. [3] are independent of N$. For 7 = 1 and F\ in the range [0.08, 0.4] the Higgs and 
the Z boson are almost degenerate in mass and so puz — 1- As a result, a m { n ~ Fi in 
this range of F\ values. Contrary to the data shown in Fig. El on Fig. 121 where a min = 1, 
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F =0.2, L = 200, y- 0.55 , on the boundary 
1 N 5 =4, M z =0.2497, M H =0.0864 ' 

n N 5 =6, M z =0.0653, M z ,=02376, M H =0.0576 
N 5 =8, M z =0.0664, M z ,=02380, M H =0.0432 
N 5 =12, M z =0.0661, M z .=0.2379, M H =0.0288 
N =16, M =0.0658, M =0.2378, M =0.0216 



F =0.2, L = 200, y= 0.55 , in the middle of the bulk 
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0.045 

^ 0.04 
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~M 
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N 5 =4, M z =0.02804, M H =0.0864 
N 5 =6, M z =0.01903, M H =0.0576 
N 5 =8, M z =0.0149, M H =0.0432 
N 5 =12, M z =0.01133, M H =0.0288 
N =16, M =0.009783, M =0.0216 




Figure 4: The combination — 0,41/' (r) + Mz / (mzr + 1) , cf. Eq. (j4.4p is plotted for different values 
of -/V5 at 7 = 0.55 and F% = 0.2 (for the boundary potential at N$ = 4 we use Mz>)- Boundary 
potential (left plot) and bulk potential (right plot). 



the values of F\ are larger than 2 and puz > 2. We compare to the results from Monte 
Carlo simulations at N 5 = 4 (diamonds) [12] and at A^ 5 = 6 (circles) [IE], using L = 12 
and T = 96. There is good agreement between the mean- field data and the Monte Carlo 
data on isotropic lattices, demonstrating that it is possible to obtain values pnz > 1- 

We have also computed the Z boson mass from the potential in the middle of the bulk 
at 7 = 1. For Fi = 0.2 (N 5 = 8), the mass Mz is the same as the one extracted from the 
boundary potential. We find pnz — 1 in the bulk for Fx = 0.2 and Fi = 1.0. 

(2) 

We find that the mass of the ground state extracted from the direct Z correlator C z 
in Eq. H3.34[) at this order in the mean-field expansion does not depend on (3 or 7 and 
most importantly it does not depend on A^. This means that using this observable one 
can measure only its infinite limit value, which for finite L, turns out to be 

Mf • = ^ (4.6) 

like on the torus [1]. This expression reflects the fact that this observable describes two 
non- interacting gluons (that is why the 4tt). 



4-3.2 Anisotropic lattice 

Motivated by Section I4TT1 we study SSB at 7 = 0.55 by extracting the Yukawa masses from 
the static potential on the boundary and in the middle of the orbifold following Eq. (14. 41) . 
As we vary N§, we set (3 to keep F\ = 0.2 constant, which means that Mh oc I/N5, cf. 
Eq. flUD. 

On the left plot in Fig. @] we show the plateaus for the boundary potential. For > 6 
there are two plateaus corresponding to masses Mz> > Mz which do not depend on N 5 . 
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orbifold, F 1= 0.2, L = 200, y= 0.55 
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mass Eq. (|4.5p in the mean-field extracted from 
) and in the bulk (right plot). 



At A^5 = 4 there is only one plateau whose value is very close to Mz> for N$ > 6. Therefore 
we identify the plateau at N 5 = 4 with the mass of a Z' boson (although we do not see the 
Z boson state). We checked that the Yukawa masses are independent of L if L is large. 
For this we compared L = 200 with L = 300 and find no differences. These data establish 
that the boundary theory is a spontaneously broken U(l) theory and by comparing to 
Eq. flOD we get a min « 0.039 N 5 for N 5 > 6. We find p HZ < 1 for N 5 > 6 and this 
particular choice of parameters as shown on the left plot of Fig. [5j Finally we remark 
that the boundary potential cannot be fitted by a four- dimensional string-like fit. 

On the right plot in Fig. H] we show the situation in the middle of the orbifold. There 
is one plateau yielding a bulk Z boson mass Mz- Mz in the bulk is decreasing as N 5 
increases and is independent on L. This result indicates that there is SSB also in the 
bulk, where we have puz > 1 as shown on the right plot of Fig. [5j The situation on 
the orbifold is therefore completely different than it is on the torus, where a Yukawa fit 
yields a mass ~ 27r/L implying that there is no SSB. This fact is by itself non trivial as it 
shows that the orbifold boundary conditions are affecting the properties of the bulk. We 
observe a difference between the Yukawa masses in the bulk as compared to those on the 
boundary. This situation is different than the one of the isotropic lattice, where we found 
the boundary and bulk Yukawa masses to be the same. 



5 Conclusions 

We have formulated the mean-field expansion in five dimensions on the lattice for an 
SU (2) gauge theory and orbifold boundary conditions. We have performed computations 
on the isotropic lattice and for anisotropy parameter 7 = 0.55 and find that the gauge 
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boson mass extracted from the static potential along four-dimensional hyperplanes is non- 
zero, both on the boundaries and in the middle of the orbifold. The gauge boson mass 
does not depend on the spatial size of the lattice, thus indicating dynamical spontaneous 
symmetry breaking. This result differs from the one obtained in the perturbative limit of 
this theory, where the gauge boson remains massless, but supports the first Monte Carlo 
simulations that were performed in [T2] . 

Acknowledgments. We thank A. Kurkela, A. Maas and P. Weisz for discussions. K. 
Y. is supported by the Marie Curie Initial Training Network STRONGnet. STRONGnet is 
funded by the European Union under Grant Agreement number 238353 (ITN STRONGnet) 
N. I. thanks the Alexander von Humboldt Foundation for support. N. I. was partially 
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A Effective mean-field action 

The effective mean-field action is defined in Eq. (12. 3p as a function of a gauge link U = 
U(n,M) and a complex matrix H = H(n,M). On the orbifold we have to distinguish 
two cases, depending whether the link is a bulk SU(2) link or a boundary U(l) link. 

A.l Bulk SU{2) links 

In order to compute Eq. (12.31) for gauge group SU(2), we start from the parametrization 

U = Uq + ia ■ u , 

H = ho-ia-h, (A.l) 

where are real and complex. It is straightforward to compute 

Retr{UH} = 2(Re h )u + (Re h A )u A 

= Reti{U(Reh ) + ia ■ (Reh)} (A.2) 

We can then define an SU{2) matrix = ((Re/i ) + icr ■ (Reh))/p with 

p = ^/(Reh y + (Reh A ) 2 , (A.3) 

in terms of which we can write 

u 2 (H)= [ DUe* B **l uir >= f Due**™. (A.4) 

JSU(2) JSU(2) 

The integral then is computed using the character expansion of the exponential 

eftrp/} = 2 _Y^{2v+l)I 2v+l {p) x ^{U) (A.5) 
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and the result is 

u 2 (H) = -h(p). (A.6) 
P 

In the mean-field we trade the 3 real degrees of freedom of SU(2) (4 minus 1, from the 
determinant constraint) for 4 independent mean-field degrees of freedom. The SU(2) 
memory is encoded in the integral above with the determinant constraint hidden in p. 



A. 2 Boundary U(l) links 

The U(l) gauge group of the boundary links is embedded in SU(2) as 

U = e^ 3 = cos(0) + i sin(0) a 3 , 

H = h -ih 3 a 3 . (A.7) 
The calculation proceeds like in the previous subsection 

Ul (H) = [ DUe kRotr{UH} = 1 d(f)e Pcos(4>) = /o(p) ) ( A _g) 

JU(1)CSU{2) 27T J 

where 

p = y/(Reh ) 2 + (Reh 3 ) 2 . (A.9) 



B The orbifold propagator 

The inverse mean-field propagator on the orbifold is 

K = - (K^y 1 + (K (vv) + K^ l) ) . (B.l) 
As on the torus, K^ vK) drops out from all expressions. The propagator is written as 

K _ x = T _ lK(hh) ^ T = -1 + K {hh \K^ vv) + K (gf) ) (B.2) 

so that it is the matrix T that is inverted when computing an observable. 

We will Fourier transform along the four dimensions but not along the fifth dimension. 
The propagator has then the components 

R- 1 = K^ip' ,n' 5 ,M' ,a';p" \ri£,M" ,a") (B.3) 

and in this Appendix we collect its pieces. 
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B.l Fourier transformation to four-dimensional momentum space 



The Fourier transformation of a double derivative of an action S 

S 2 S 



Sv'Sv" 



with respect to v' = v a i(n', M'), v" = Va»{n", M") defines a kernel 

K(p',n's,M',a'- P \nlM",a"). 
The Fourier transformation is 



r (4) 



5 2 S 



(B.4) 



(B.5) 



(B.6) 



[ %) W 5v'5v"' 

where = TL 3 is the four- dimensional volume. The factors (— i) s M'& and (i) 5 "" 5 are 
introduced to make the propagator real, a property that reflects CP invariance [23J. 

The kernel Eq. ( IB. 51) is a matrix which is divided into sub-matrices in coordinate 
indices (n' 5 , n'l) for given (p', M', a'; p", M", a"). These submatrices have dimension (N$ + 
1) x (JV fi + 1) if M' = // and M" = //"; dimension iV 5 x (JV 5 + 1) if M' = 5 and M" = /i"; 
dimension iV 5 x iV 5 if M' = 5 and M" = 5. 



B.2 f^C 1 ^ in four-dimensional momentum space 

on the orbifold is defined as the second derivative 

a 2 



dh a >(ri,M')dh a »(n",M") 
of the orbifold effective action. We introduce the notation 



a 2 (a, 0) = a 2 (a, iV 5 ) 



if a = 1,2 
a 2 if a = 3 



and use the quantities 

Oi 

MO 



llu'M (u' x {p)\ 
2 I P P \ P 



p=h (0)=h (N 5 ) 



2 1 p 



p=fto(0)=/ l o(Afs) 



^ 2 (p) ! / ?4(p)V 

P VP 
« 2 (p) 



p=M0 



(B.7) 



(B.8) 



P 



(B.9) 

(B.10) 

(B.H) 
(B.12) 



P=ho(-) 
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As mentioned, we perform a Fourier transformation along only the four-dimensional hy- 
perplanes. The a' = a" = component is 

{JVb-1 n 6 -i 
(1 - &M'b) ^ S n' 5 ,m 5 bl(m 5 ) + S M '5 ^ ^,m 5 &l(^5 + V 2 ) 
ms=l m,5=0 

(B.13) 



and 



K (hh),bound. = S^ p „5 M , M „5 n >> jn > 5 (l - 8 M 'b) {8n' B ,0 + ?>n' B ,N 6 ) «1 ( B -14) 



The a' = a" = A ^ components are 



K (hh),bulk = 5^ pll 5 M , M „5 n ,, nl ) (1 - 5 M , 5 ) 5 n ^ m5 b 2 (m 5 ) + 5 M '5 ^2 s n' 5 ,m 5 b 2 ( m 5 + V 2 ) 



JVs-l 



AT 5 -1 



m 5 =l 



m,5=0 



(B.15) 



and 



^ftj),bound. = 5^ pll 5 M , M „5 n ,, ^(i - S M , 5 ) (5 n > 5i0 a 2 (A, 0) + 8 n ^ Nz a 2 (A, N 5 )) (B.16) 



B.3 in four-dimensional momentum space 

B. 3. 1 Gauge fixing 

We use the backward derivatives f A (n, M) = v A (n, M) - v A (n - M, M), A= 1, 2, 3 and 
introduce the gauge fixing term 

2 



' N 5 -l 



+ z l 



{n M } I. «5=1 A 



E/ A (n,/i)+ 7 />,5) 
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^/ 3 (K}^5 = Ar 5 ^) 



We set the boundary weight to 



1 

Zl= 2 



The Fourier transformation of (v' = va'(ti', M'), t>" = VA"(n", M")) 



we denote it by 



5v'5v" 

K^(p',n' 5 ,M',A';p",nlM",A") 



(B.17) 

(B.18) 

(B.19) 
(B.20) 
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and is divided into three contributions. 
Contribution 1 is for W = //, M" = //': 



, N 5 -l 
^ I A n 5 =l 



+ z 1 5 A '' A "5 A ''X' 5 ,o + zj**" 6**6^ , (B.21) 

where ]3 M = 2sin(p M /2). 

Contribution 2 is for M' — 5, M" — //': 

_ N B -1 

-r$W»" SA '' A " E ^ E l 6 n>MK - <Wi<Wi] ■ (B.22) 

? A n 5 =l 

The contribution for M' = //, M" = 5 is obtained using the Hermiticity property. 
Contribution 3 is for M' = 5, M" = 5: 

(B.23) 

^ A n 5 =l 

where 

□n^X 5 '( n 5) = S n > 5>n »(5 n ' 5>n5 + 5n^n 5 -l) ~ ^<,< + l^,n 5 -l - ^'.r^-l^.ns ■ (B.24) 

We have also implemented a background gauge fixing 

f A (n,fi) = v (n 5 ) [v A (n,fi) - v A (n - fi, fi)] , 

/ A (n,5) = v (n 5 -l)v A (n,5)-v (n 5 + ^)v A (n-5,5), (B.25) 

which is equally valid and vanishes when evaluated on the background v . It does not 
change any of the results for physical observables. 

B.3.2 The Faddeev-Popov determinant 

Even though not directly relevant for the gauge propagator we can now carry out the 
Faddeev-Popov construction, necessary for the free energy. The ghost action is 

fe=EE c A VMtU" cA V) (B.26) 

A', A" n',n" 

where the sums run in the fundamental domain of the orbifold. The ghost kernel 

M A'n';A"n" = E M A^n';A"n" (B.27) 
M 
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is obtained from the variation under infinitesimal gauge transformations 

8f A \n,M)= Y, M *n-,A>>n^ A '\n"). 



(B.28) 



A",n" 



The changes on the orbifold with respect to the torus case in [I] is that on the boundaries 
the gauge transformation is U(l) and acts only on v 3 at the boundaries at = 0, N$. 

We define the matrices 



8 ' ^0(^5) (25 n ",n' — 5n",n'+u ~ ^n",n'-v) 



with four-dimensional Fourier transformation 



M {v) 



(B.29) 



;B.30) 



and 



M (5) 

JVl A'n' b ;A"n'l 



5 A ' A " ( v (n' 5 + h+ - 1) ) £„« 



S A ' ,A " + S^.J*' A "(5 A '> 3 - l)j v (n 5 + ^) %,^ + § 



5 



A' A" 



1 



(B.31) 



with the four- dimensional Fourier transformation changing <L« „/ — > 5^ „■ The final 
expression for the Faddeev-Popov kernel is 



E 

Z 5 =l 



+ 



*1 E [^.O-Mg-S.n" + 5 <,^-^3X; 3/ , 



(B.32) 



The relevant for the free energy Faddeev-Popov determinant is 

D ( °; h) = detM {orh) . 



(B.33) 



B.3.3 The double derivative of the plaquette action 

As for K^ hh ' due to the gauge invariance of the Wilson plaquette action the contributions 
of double derivatives with respect to links to the orbifold propagator simplifies to 



dv a >(ri,M')dv a/ ,(n",M") 
+ boundary terms 



~H~T7ii\ E E Sw[U(n, M] 



neT M 



(B.34) 
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with F the fundamental domain of the orbifold with the boundary contributions sepa- 
rated out. 

The gauge structure of the full kernel is 

° 1 (B35) 
Ka ' a "~ { (A',A") 

where we use the notation (A', A") for the blocks along the algebra indices. Since 

S 2 S W 5*S W (B36) 



Sv' Sv" Sv" Sv' ' 

the kernel is Hermitian. Hence, the off-diagonal elements are related by the Her- 

miticity relations 

feO ={{K%1) V (B.37) 

V J M'M" \\ J M"M'J 

where ( K^, v 2,) is a (N5 + l)-dimensional matrix. 

V / M'M" 

j((vv) Jepenjg on the following weights 

= 5«', Q "^o(^5 + l/2)v (n 5 + 1) , 
: S a ',a" v (n 5 )v (n 5 + 1/2) , 

(a',a")/ \ 



(a' 


,a"), 




(a' 


,a") ( 


fa) 


(a' 


,«"), 




(a' 
«2 


.a"), 





v (n 5 - l/2)v (n 5 - 1) . (B.38) 



and 



w[ a,a \n 5 ) = S a >, a „(v Q {n b )) 2 , 

w ( 2 a ' a \n 5 ) = S a/ ^v (n 5 )v fa + 1) , 

w ( 3 a ' a \n 5 ) = 5 a > >a » {v fa + 1/2)) 2 , 

w ( f' a "\n 5 ) = 5 a ,, a „{v Q {n b -l/2)) 2 . (B.39) 

The background is defined only for Vofa), ^5 = 0, 1/2, 1, 3/2, • • • , N 5 . Outside this range 
it vanishes by definition. Correspondingly, the above definitions of the weights hold for 
every value of 715 = 0, • • • , N 5 in the range where the background is defined and for all 
values of (a', a") except from a few special cases which are related to certain contributions 
from boundary or bulk/boundary plaquettes. We specify these special cases below. 

To begin, we set for A' = A" = 1 and A' = A" = 2 the weights that originate by 
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taking at least one derivative on one of the boundaries to zero: 

^(O) = w?' 1) (N 5 ) = 0, 
w^\N 5 ) = 0, 

«? ,1} (o) = o, 

^) (i V 5 _i) = o, 

u^\0) = u£> 1 \N 5 )=0, 

w^\l) = ^ 1 ' 1) (iV 5 -l) = 0, 

w { 3 ,l) (0) = 0. (B.40) 
For the elements A' = A" = 3 we have 

w (3.3) (0) = 1/2 ^( ) 2 , 

^f' 3) (A^ 5 ) = ^o(iV 5 ) 2 , 

w ( l' 3) (N 5 ) = v (N 5 -l/2)\ 

4 3,3) (iV 5 ) = v (N 5 -l/2)v (N 5 -l), 

W ^\N 5 -1) = v (N 5 -l/2)\ 

4 3 ' 3) (iV 5 -l) = v (N 5 - l)v {N 5 - 1/2) . (B.41) 

As mentioned, all other weights have their usual value, defined in Eq. flB.38j) and Eq. (IB.39j) . 
Weights that are not defined are assumed to be identically zero. Finally, the anisotropy 
can be also absorbed in a redefinition of the weights: 

(a',a") 
V l,2 ->■ u l,27> 

(o',o") 



(a', a") 
1 



— > ifi/7 



^w' ->• ^2,3,47 (B.42) 

for every (a', a"). 

We will now compute the various individual contributions to the full kernel. We define 
the symbols 

e(a'; M\ M") = sgn [W» + (1 - <W0(2<^,o - 1)] (B.43) 

and 

e(p',n' 5 ,Af',a';p",n' 5 ',M",o")(^ / ) -W' > a i C i <$) = 
t ii 

(a8 <t< + e(a'; M', M") &<^+v 5 e ~^ ,(1 ~ % ' s) 
+ e{a' ] M\M'')c5 < ^ 5N , l ^^ 

(BAA) 
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The double derivatives contribute then the components 



K^ vv \x'; x") = -0. 



Sm'm"Sm'5 e( x / ;a ,//) («, fi] 0, w 2 , w 2 , 0) 



+ 5m'M"(1 -S m > 5 )( ^2 e {x';x")(fi>, MOiW^w^O) +e( a ./ ;x »)(5,5;0,w 3 ,w 4 ,0 



+ (1 - 5m'M")(! - £m'5)(1 - ^iw»5)e( a: '; a; '/)(M / ,M";wi,w 1 ,wi,wi) 
+ <W 5 (1 - ^Af"5)e(x' ; z") (5, M"; v u v 2 , v 1 ,v 2 ) 

+ (1 - 5 M r 5 )8 M "5e(x';x") { M> \ 5 ; Ml, Ml, U 2 , U 2 ) 



(B.45) 



The argument of all weights in the above is n' 5 and their superscript (a', a"). Also we 
have defined the collective index (V; x") = (j/, n' 5 , M', a'; p", n' 5 ', M", a"). 

To these components, the corresponding components of the contributions from the 
gauge fixing term must be added, so that the full contribution from the pure gauge part 
of the action to the propagator is 

K {vv) + iT (gf) . (B.46) 



We now have all the ingredients to compute K 1 in Eq. flB.2j) . The last issue to be 
discussed before one does so is its eigenvalue structure. 

A general property of T is that its eigenvalues are invariant under — > —p^ for any 
fi = 0,1, 2, 3. This is useful when computing the free energy numerically. Zero eigenvalues 
in K^ hh \ if any, clearly do not contribute anything to the free energy or to any of the other 
observables. The matrix T may have certain zero eigenvalues which on the other hand 
must be taken care of otherwise it cannot be inverted. We separate these zero modes in 
two classes. 

One obtains spurious zero eigenvalues due to the Dirichlet boundary conditions that 
force some of the fields to vanish on the boundaries. For example, for a' = a" = 1 and 
a' = a" = 2 the link variables v' a (0, fi), v' a (N 5 , fi) all vanish resulting into 2(d — 1) (d = 5) 
zero eigenvalues in the 11 and 22 components of T. These spurious eigenvalues occur for 
every value of the four-dimensional momenta and appear just because bulk and boundary 
components are packed together into the propagator. They are unphysical and can be 
removed by hand. 

The second class of zero eigenvalues of T contains those corresponding to vanishing 
four-dimensional momentum p = 0. On a finite lattice with fully periodic boundary 
condition they appear in any gauge invariant formulation and persist even after fixing 
completely the gauge, due to a left over global gauge invariance surviving on a lattice [2]. 
Beyond the spurious eigenvalues mentioned above, the following properties hold for any 
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a 7^ if = V fi: without gauge fixing T has N 5 zero eigenvalues, corresponding to local 
gauge transformations of the links along the extra dimension. After gauge fixing through 
Eq. ( IB. 171) (also in the variant of Eq. (IB. 251) ). only one zero mode survives. This is called 
"toron" . If p ytz then there are N5 — 1 zero eigenvalue which are completely removed by 
the gauge fixing. When present, such toron zero modes render the corresponding physical 
observable, plagued by infinities, unusable. Since however it is a finite volume effect, if 
(and only if) one can ensure the existence of some regularization which leaves behind a 
finite result, they can be dropped; their contribution to any physical observable is volume 
suppressed and disappears when the infinite volume is taken pp. 

We finally report on the zero eigenvalues of the Faddeev-Popov determinant. When 
the four dimensional momentum p vanishes there are 2 zero eigenvalues for any a. They 
correspond to lines of zeros in the Faddeev-Popov matrix Ai. If p is not zero there are 
two zero eigenvalues for a = 1, 2 and none for a = 3. 



C Z and Higgs to leading order 



We first define quantities that will be heavily used in the calculation of the observables. 
Let 

A^K) = ^ -^— = (1-6^) — 1 (C.l) 
f^ Q u„(r + l/2) Va(n 5 + ±) 

and 



0) = 2^cosj9' t Ai iY5) K)^ M) (Po,0,n' 5 ,5,a;p' ,0,n / 5 ',5,/3)AS 



PO n 5' n 5 



n 5, 



PO n 5> n 5 



(C.2) 



We start from the Z mass correlator. The computation is similar to the one for the torus, 
so for more details see pp. First we look at the single derivative of — p(°)t which after 
the derivatives leaves behind the matrices 

T? = a a gg* + go a gi-h.c., (C.3) 

in the background. A simple calculation gives 

S° = 0, S 1 = Aia\ S 2 = Aia 2 , S 3 = . (C.4) 

The group structure contributes terms of the form tr{a 3 E Qi ^E Qj }, to be contracted against 
the Euclidean structure, which we show below: 
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times 

eW-ptis^S^S^S^K-^at!, a 3 )K~\a 2 , a 4 ) 

+ e Wi+P' k )§^ 25 a 2 A l5 a 3 A 35 a 4 A 4K -l^ a ^ R ^^ a J 

+ e -i(p|+p' fe )^i^i^A 2 r A3 r i4 r i| ai) a 3 )ir- 1 (a 2 , a 4 ) 
+ e-*^-^^ 01 ^^ 02 ^ 1 ^ 03 ^ 4 ^" 4 ^^-^^!, a 3 )K-\a 2 , a 4 ) 

+ e ^' i < i 6 aiAl 6 oaAa 6 aaAii 6 a * A *K- 1 (a 1 , a 4 )K-\a 2 , a 3 ) 

+ e-'W-K-)^ ^ 2 ^^ 03 ^^ 4 ^- 1 ^!,^)^" 1 ^,^)- (C.6) 

Putting everything together, the observable can be first written as 

Cf{t) = ^(^o (0) ) 4 (-o(0)) 4 -^ £ tr{a 3 S^t S A 2}tr{ -3 S A 3ts A 4} . 2 . j 

^ ' Ai,-,A 4 p" k,l I 

+ e^+^n^^, As)^^, Ai) + e-^^ng^CAt, A 4 )ng 1) (A 2 , A3) J (C.7) 
and finally be simplified to 

= £ Sr^ ( p o 0) ) 4 (^(0)) 4 £ E ■ nj;; i)( i, i)ng„(2, 2). (c.8) 

We denote by Pq ^ the Polyakov loop evaluated on the background. In the above we have 
used that the elements of the propagator which are diagonal in the M index are invariant 
under p 1 — > —p 1 . The coefficient is the symmetry factor of the double exchange 
diagram and it is 

<* = ( 4 ) . (C,) 

We note that the toron (p^ = Wfi) does not contribute to Eq. flC.8j) . 
Next we turn to the Higgs mass. The relevant trace is now 

{i)a = to{o a gg* + go a g*}. (C.io) 

Apart from the traces originating from the above structure, the computation is again 
analogous to the one for the torus geometry (see p^) so we give directly the result: 

c&(t) = ^(^o (0) ) 2 ic (1) E< 1 >- II S 1 >( a ' a )' ( c - n ) 

a 
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summation over repeated gauge indices implied. The symmetry factor is 

c« = l. (C.12) 

Performing the traces one obtains 

^ ) w = ^iy(^o (0) ) 2 ^n;;; 1) (o,o). (c.i 3 ) 
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